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Abstract 

System D 2 y(t) + z(t) = f(t,y), f Q 2 M~f)D"<y(t)dj = f Q 2 M~t)D"< z(t)dj, t > 0, where D 7 , 7 £ 
[0, 2] are fractional derivatives, is investigated and the existence of the solution in a mild and 
classical sense is proven. Such a system arises as a distributed derivative model of viscoelastic 
body and in the system identification theory. Also, the condition for existence and unique- 
ness of a solution to a general linear fractional differential equation Y^i=i o-iD^zit) = g(t), 
Yi € M, i = 1, k in y+ is given. 
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1 Introduction 

In this paper, we solve and analyze solutions to a nonlinear system motivated by a mechanical 
model 

D 2 y(t) + z(t) = f(t,y(t)), t>0, 

< r ' 7 '> i fa MWviW-y = Jo M7)^z(t)d 7 , t>0, 

y(0)=yo; y ( - 1 H0) = v Q . 

Here, y and z, which represent strain and stress respectively, are unknown functions defined 
for t > 0. Given are locally integrable function / representing forcing term, and (f>i,(f>2 which 
characterize the material under consideration. Constants ya and "Uo are initial displacement and 
velocity. Operators D 1 , 7 e [0, 2] are operators of usual (for 7 integer) or fractional differentiation. 
System (CP) models various physical processes. For example, (CP)\ represents an equation of 
motion of a physical pendulum (in this case f(t, y(t)) = siny(t)), moving in a dissipative media of 
viscoelastic type; (CP) 2 represents a constitutive equation of a material where energy is dissipated. 
Since the upper bound in integrals in (CP)2 is two, both viscoelastic and viscoinertial effects are 
included. As experiments show, the presence of integral on the left hand side of (CP)2 indicates 
that dissipation properties depend on the order of the derivative. The integral on the right hand 
side of (CP) 2 is a consequence of the well known principle of equipresence (cf. 23j ) . 
Fractional operators D 1 and / 7 , 7 g I, (see ^ and (H), respectively) are widely used in consti- 
tutive equations for modeling viscoelastic bodies. A typical model is the generalized Zener model, 
where relations between stress and strain involve linear fractional differential operators. We refer 
to [2], [H] and [8] for details. A viscoelastic model of wave propagation with distributed order 
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derivatives was analyzed in [Uj. In [T] a distributed derivative model of viscoelastic body was 
proposed, and the restrictions which follow from the Second Law of Thermodynamics are stated. 
Distributed order derivatives appear in many other branches of physics. For example, Caputo in 
10J introduced distributed order fractional derivative to model dielectric induction and diffusion. 
In the series of papers [HI [T31 E3 HO] the fractional diffusion equation, which was obtained from 
standard diffusion equations by replacing the first-order time derivative with the fractional deriva- 
tive of order £ (0, 1) was analyzed. In [7] and [15] distributed derivative models were used in 
the context of system identification theory. 

Nahusev, [21] , considers distributed order derivatives of the form D a yda, c,d£l. He used the 
Laplace transform method to show the existence of a solution to equations of type f D a yda = g, 

c < and /„ D a yda = g, d > in L l [a, b], with g S L l . Recently, Kochubei [H] gave the 
necessary condition on continuous function <fi in order that distributed order fractional differential 
equation of the form J Q (f>(j)D a yda = g has a solution. Some examples of mechanical models 
with distributed order fractional derivative were also studied in our previous works [3] , [3] and [5] . 
Those are all particular cases of system (CP). 

In order to solve system (CP) we will use the Laplace transform of tempered distributions sup- 
ported by the half line and the Schauder fixed point theorem which leads to a solution in an 
interval [0, S], for some 8 > 0. Also, we will show existence and uniqueness for a distributed order 
fractional differential equation J , (/>(7)Z? 7 z = g in if <j> is of the form <j> = ^2 i=0 ai§(- — 7i), 
a, G R, di 0, 7i G R, i = 0, 1, k so that this equation has a form of a linear fractional dif- 
ferential equation ^i=o a iD li z = g. We refer to [T7] for solutions of linear fractional differential 
equation in spaces of continuous or locally integrable functions. Also we refer to |16j for solutions 
within spaces of distributions defined on the hole line obtained by the use of the Fourier transform 
method. 

In Section [21 we will define different types of solutions to (CP) problem: classical, mild and 
non-impact ones. Further on, in Section O we state and prove the existence and uniqueness of a 
solution to a general linear fractional differential equation in the frame of , assuming that the 

corresponding holomorphic function X^i=o ais7i ' < 7i < 2, Res > 0, does not have zeros. The 
existence of solutions to system (CP) is proved in Section 2] In Theorem 14. II we give sufficient 
conditions for the existence of mild, non-impact and classical solution to (CP), in the case when 
4>i is continuous in [c, d] for some < c < d < 2 and ^(O = Y^i=o a «^(' — 7*)j < 7i < 2. Further 
on, we analyze cases where form of <f>2 remains the same, while <j>\ changes. In Section [5j we 
perform our analysis for the case ^1(7) = b 1 and ^2(7) = a7 , b > a, where b > a is a restriction 
which follows from the Second Law of Thermodynamics. 

1.1 Notation and notions 

Let I C K be an interval of K. Then, we denote the space of locally integrable functions in / by 
Lj (J), the space of fc-times continuously diffcrcntiablc functions y by C k (I), and the space of 
functions y G C^ 1 ^), such that G L\ oc (I) by AC k (I). 

The space of functions y G Lj oc (E), with the property y(t) = 0, t < we will denote by L 1 1 oc+ (IR). 
Further, £^(R) and o5^(R) are the spaces of compactly supported smooth functions and of rapidly 
decreasing smooth functions in R, respectively. Their duals 5^'(R) and ^'(R) are the spaces of 
Schwartz's and tempered distributions, respectively. f^_(R) and J^f(R) denote their subspaces 
consisting of distributions supported by [0, 00). Henceforth, we drop R in the notation of @' + (M) 
and c5^(R). We will also use compactly supported distributions, elements of £'(ft), where fl is 
an open interval of R. Any element of £'(f2) can be extended to R as an element of £^'(R) in the 
usual way. Let J C I be a subinterval of / and let y be a function or a distribution defined in I. 
We will use the same notation y for its restriction in J. Thus y\j will be denoted by y again. In 
this sense, we will write / G S?'+ n C p ([0, a}), which means that / e and /|[ , o ] G C p ([0, a}). 
The Riemann-Liouville operator of fractional differentiation or integration is denoted by D 7 , 7 G R 
defined as follows. Let y G L 1 1 oc+ (R). The Riemann-Liouville fractional integral of order 7 > is 
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defined by 



where T is the Euler gamma function. For 7 = one defines I°y := y. If y £ L 1 1 oc+ (R), it is shown 
(cf. [22l Theorem 2.6]) that lim^^o Hy(t) = y(t) almost everywhere in R. 

Let y £ L^^R) and suppose that for every a > 0, y £ AC k ([0, a]). The Riemann-Liouville 
fractional derivative of order < 7 < k for some k £ N, is defined by 

Jk 

Wy(t) := -^I k ^y(t), t > 0. (2) 

If 7 £ N, then £H = Jl. 

Ricmann-Liouville's fractional operators in the setting of distributions are defined (e.g. [23]) by 
introducing the family f a £ , «el 



/a(i) := 



ff(*)l75j» i£R, a>0, 



where iJ is Heaviside's function. Then, f a * fp — f a +p, for all a,/?£K and f a * is the convolution 
operator in Also, f a : Sf[ — > tP^.. The operator / Q * in is the operator of fractional 

differentiation for a < and of fractional integration for a > 0. It coincides with the operator of 
derivation for —a £ N and integration for a £ No = N U {0}. In particular, /o = S. 
If a > and y 6 Aoc+W then I a y = f a * y. If < a < fc, k £ N and for all a > 0, 
y £ L 1 1 oc+ (M) n AC k ([0, a]), then Z) Q y = /_ Q * y. We will use notion D a y also when y is distribu- 
tion and will denote D~ 1 y = 7 7 y, 7 > 0. Note that D 7 / 7 y = y for y £ Lj^^R) and l 1 D 1 y = y, 
7 > if y is considered as an element in &' + . 

Let < c < d, 4> e C([c, d]) and y £ L 1 1 oc+ (R) n AC 2 ([0, a]). Then the distributed order fractional 
derivative in [c, d] is given by 

<K7)i^y(i)d7- (3) 

If </> : [c, d] — ► R+ be a piecewise continuous bounded function and y £ J5^_f_ , then as in [5] , one can 
consider ^ as a Bochner type of integral. We refer to [5] for the properties of © within 
Recall the definition from [5]: 

Definition 1.1. Let £ £'(R) and y £ ^_|_. Then J (p(^)D 1 y is defined as an element of 

supp 

by 

' 4>{i)jyy{t) dry, <p(t)) := (0( 7 ), (^l/W, ¥>(*)», V £ ^W- 



When supp</> C [a,b] we write J b 4>{"i) D 1 y dj instead of J </>(7) D 1 y dj . It is shown in [5] that 

supp <fi 

7 1— > (£> 7 y, 4>) : R — > R is smooth and that the mapping (a, y) 1— > D a y is continuous from R x J^l. 
to =5^. Also, (/ ^ J <j)(~f)D 1 y d^ is a linear and continuous mapping from ,9"^ to S?\. For 

supp (p 

k 

M-) = -7*). 7i€Ri£{0,l,..,fc} (4) 

z=0 

distributed order fractional derivatives becomes a sum of finite number derivatives of fractional 
order £\ =0 OjD^y. 
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Distributed order fractional differential equation is given by 

J ct>( 1 )D^yd 1 = g, (5) 

supp <f> 

while for tj> of the form it becomes linear fractional differential equation. 

To deal with fractional differential equations, it is convenient to use the Laplace transform. Recall, 

if y G y! ll _ 1 then its Laplace transform is defined by 

y(s) = Cy(s) = (y(t), ^(t) e - st ), Re a > 0, 

where tp G C°°, tp = 1 in (—a, oo) and tp = in (— oo, —2a), a > 0. Note that £y is an analytic 
function for Re s > 0, and that the definition of Cy does not depend on a chosen function tp with 
given properties. We will use the identities given in the next proposition (see [5]). 

Proposition 1.2. Let 4> G £'(M) and y G Then: 

a) C(f a * y)(s) = ±y(s), Res > 0, a G R. 

b) C( J ^)D~'yd 1 )(s)^y(s)( ( j>( 1 ),s''),Res>0. 

supp 

c) If (j) is a continuous function in [c, d] and ^(7) = 0, 7 ^ [c, d] ; then 

£( 0(7)i^»d7)(a)=*(*) / 0(7)s 7 d7, Re.s>0. 



2 Definitions of solutions to (CP) 

Problem (CP) will be analyzed in the context of the following definitions. 

Definition 2.1 (Classical solution). Let f : [0, oo) x R — > R 6e continuous and fa, fa G £'(1 
supp& C [0,2], i = 0,1. Lety ,v a G R, <5 > 0, y G y|nC 2 ([0,^), z e y{nC([0,(S]). 
We say i/iai (y,z) is a classical solution to (CP) in [0,5] if y and z satisfy (CP), where 



(CP) 



D 2 y(t) + z(t) = f(t,y(t)), te[0,6], 
J 2 fa( 1 )D''yd 1 = J 2 faWHrzdry, in 
2/(0) =2/o, tfW(0)=«o. 



Let fa, fa G £'(R), assume that y G ^+ and let z G J^l. be a solution to ([5]) with (j) ~ fa and 
.9 = J Q 2 4>i(-/)D' 1 yd'y. Then z is a solution to (CP) 2 in This is provided by 

1 



(62(7), *» 7 > 



and then 



1 



(6) 
(7) 



,(02(7),s 7 ), 

Equation (CP)\ with initial conditions (CP)3, in the setting of tempered distributions, reads 

/_2*tf + * = /(-,i0 + i/b<5' + «o& (8) 
We substitute in © and apply D~ 2 = f 2 *. ft follows that 

* / 1 ( 7 ) J D 7 - 2 yd 7 + i?- 2 /(-,y)) + «o/i + ^. 



(fa (7), « 7 ) 
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Definition 2.2 (Mild solution). Let f : [0, oo) x K — * R be a locally integrable function and 
4>i,4>2 S £'(R), supp fa C [0, 2], i = 0, 1, where 4>2 is such that (0|) holds. Let yo,Vo el emd (5 > 0. 
The pair (y,z) is called a mild solution to (CP) in [0,5] if y S R C([0, <5]), smc/i i/iaf 

+ ^- 2 /(t,?/(t)) + ¥ + j;o ! te[0,i], (9) 

and z is given by Qj. 

Remark 2.3. a) A mild solution to (CP) in [0,(5] is a weak solution in (0,6). Namely, © and 
© imply that for every 9 G C °°((0, (5)) 

(£> 2 y, 0) = (y, D 2 9) = (D- 2 z(t) + D~ 2 f(t, y(t)) + v Q t + y , D 2 d(t)) 
= (z,9) + (f(t,y(t)),9(t)). 

b) A classical solution is also a mild solution. If / : [0, oo) x R — > R is continuous and a mild 
solution exists, then considering that mild solution satisfies (CP)2, in order for it to be a classical 
solution, one needs additionally to assume that y € C 2 ([0, 6]) and z € C([0, <5]). 

In mechanics, an impact solution is an absolutely continuous function whose first derivative has 
jumps. Impact solutions are included in the definition of the mild solution. In order to make 
an additional classification of solutions related to the observable quantity y, we introduce a class 
of non-impact solutions. This class will correspond to solutions in which the first derivative y' 
(velocity) is continuous. 

Definition 2.4 (Non-impact solution). Let f : [0, oo) xl^l be a locally integrable function, 
<^i,<^2 £ £'(K), supp</>i C [0,2], z = 0,1. Let yo,vo £ R and 5 > 0. The pair (y,z) is called a 
non-impact solution to (CP) in [0,6] if y e n AC 2 ([0,6]) and z 6 ,9" + V\L x ([0,S\), such 
that 

D 2 y(t) + z(t) — f(t,y(t)), almost everywhere t G [0,6], 

Ml)D^y(t)dj= f 2 ( 7 )^z(t)d 7 , m y\_ 
o Jo 

and 

2/(0) = y , Dy(0) = v . 

It is clear that condition ([6]) plays an important role in definitions of the solutions, and that 
satisfaction of © depends on <\>2- However, in the next section, we will see that when </>2 is of the 
form © with additional condition (An) satisfied, then © is fulfilled. As we will sec in Section [H 
© and certain conditions on / assure the existence of mild, non-impact, and classical solutions. 



3 Linear fractional differential equation in 5?\_ 

In this section, we prove the existence and the uniqueness of a fundamental solution to a general 
linear fractional differential equation in S^[. 

Let g G Let ^ and 7, be arbitrary real numbers, such that 74 > 74+1, i G {0, l...,k}, and 
consider 

k 

^a.D^z = g, in Y[. (10) 

i=0 

Theorem 3.1. Equation M0\) has a unique solution z G 5?L , if and only if 

k 

(A ) ^a^^O, s G C+ = {s G C; Res > 0}. 

i=0 
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Proof. First note 

y^gjD^z = ^ai(/_ 7i *z) = (^a»/_ 74 ) * z, 

t=0 i=0 i=0 

i.e. that (|10[) is equation of a convolution type. For given ip, h G J^f, the equation ^ * 2 = h is 
uniquely solvable if and only if there exists a, b G R and C > 0, such that 

i^- ,, < C— r- r7—, s G (L+. 

|^(s)| - |Res| & 

where ^ denotes Laplace transform of ip, (see [25, Chap. 7, p. 50]). 
Let us denote 

fe 

F(s) = 0(s) = ^a lS 7 \ s G C+. 

i=0 

We will show that F(s) ^ 0, s 6 C+ (which is (Aq)) implies that there exist a€l and C > 0, 
such that 

1 (l + |s|) a 

This will imply the theorem. 

To begin with, note that function F is analytic in C\(— oo, 0]. Next observe that | ^(s) | behaves 
like aos 70 1 when |s| is large, and like \a,kS lk | when s is small which implies that there are r, R G R, 
< r < i?, such that F(s) has no zeros for \s\ < r and \s\ > R. Further observe that F(s) ^ 
for Res > implies that F(s) ^ for Res < and therefore the set of zeros of function F is 
contained in [—iR,—ir] U [ir, iR]. Thus, it must be a finite set of points iyj, j G {1,2. .J}, since 
zero set of an analytic function does not have limits in the domain of its analyticity, so in any 
compact region, the set of zeros can only be a finite set. 

Since F is analytic, it follows that in neighborhood of zeros iyj, \F(s)\ behaves as |s — iyj\ mj , 
rrij G N, j G {1,2, ..,/}. Denote Dj = {s G C; |s - iyA < rj, for some r, > 0, j G {1,2. ..I}, 
D = U^ =1 Dj and K = {s G C;r < |s| < i?}. Choose a > rrij, j G {1,2...,/} and a > 7o . There 
exist C such that 

— \—r < C , 1 r , s G (C\if) U £>. 
|F(s)| " (Res) a ' V X ; 

In K\D function F reaches its minimum, denoted by m, so with some new C 

— — / 1 ^ ... ' 



\F(s)\ ~ m ~ (Res) a 
Putting all together we arrive to (fTTj) . □ 

Remark 3.2. Equation (JTQJ) represents a constitutive equation of a viscoelastic body which, 
coupled with equation (CP)\, represents a mechanical model (CP). We explained in [S] the 
mechanical aspect of condition (Ao) and here we repeat it briefly. Let T > and y be of the 
form y(t) = J2i=a a i D li z(t), t G [0,T], The dissipation work is given by Ad — J Q z(t)y^ (t)dt 
and the dissipation inequality reads Ad > (see [S|). If (Ao) does not hold, then the dissipation 
inequality is violated. So, (Aq) is a necessary condition if (CP) represents a mechanical model. 

The solution to (fTI)]) is of the form z — l*g, where I is a fundamental solution to (|10p . i.e. solution 
to y\._ aiD li y = 8, with the properties given in next lemma. For the proof see [5]. 

Lemma 3.3. Let 7.; G [0, 2) and 7, > 7 j_|_i > /or all i G {0, 1..., fc}. Assume (Aq). Let 



!(*)= £ te^J W ' ^° (12) 
, t < 0. 



TTien: 

(i) I is a locally integrable function in R. 

(«J Z is absolutely continuous in R, j/70 — 7 fc > 1. 
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4 Solvability of the system (CP) 

Through this section we consider cases when 02 is of the following form 

fc 

M-) =^2a i 8(---y i ), 7 i e[0,2),7 < >7 i+1 ,te {0,1, ...,*}. (13) 

i=0 

The constitutive equation (CP) 2 then reads 

r-2 k 

/ 4> 1 ( 1 )D~'y(t)d 1 = Y^a l D^z(t), t > 0. (14) 

If condition (Aq) of Theorem l3.1l is satisfied, and I is defined as in Lemma |3~31 by (H2), then, if we 
suppose y G J^l. , Theorem 13.11 implies that (fl4"|) has a unique solution z G =5^ given by 

/• 2 

z = Z * / tfii(j)D' r yd'-f. 
Jo 

Equation @ becomes 

y(t) = -l* f Ml)D"<- 2 y{t)dl + D- 2 f(t,y{t))+v Q t + y Q , t > 0, (15) 
Jo 

so a mild solution to (CP) will exist if (|15p has a solution t/ G C([0, $]) fl =5^, for some <5 > 0. To 
show such existence, we will define a mapping T by the right hand side of (|15|i , and show existence 
of its fixed point in appropriate Banach space. The necessary conditions for T to have a fix point 
are to follow. 

First, let us impose conditions on /: 

(Ai) f(t,u), (t,u) G [0,oo) x R, is locally integrable and 

there exist a > 0, a > and h G L ([0, a]), h> 0, such that 
|/(t,u)| < M*)M a , u G R,t G [0,a], 

and the stronger one 

(Ai)' f(t,u), (t,u) G [0,oo) x R, is continuous f(t,0) = 0, £ G [0, 00), and 

there exist a > 0, a > 0, h G C([0, a]) with h>0 and /i(0) = 0, such that 

\f(t,u)- f(t,v)\ <h(t)\u-v\, u,veR,te [0,0]. 

Further, we introduce conditions on 0i and 02- 

($1) 0i is continuous in [c, d] C [0, 2), ^1(7) = 0, 7 ^ [c, d]; 
</>2 is of the form (fTSl . 

The following theorem holds. 

Theorem 4.1. Lei j/o,«o G R, r > max{|?/o|, \vo\}- Let ($1), (Ao) and (Ai) hold. There exists 
6 = S(r) > such that: 

a) (CP) has a mild solution in [0,6]. 

b) If y(0) = 2/0 = 0, i/ien (CP) has a mild solution (y,z) in [0,5], such that y G C^QO, 1]). // 
7o — Ik > 1 ifeen (2/, z) is a non-impact solution. 

c) If (Ai)' holds and yo = Vq = 0, then (CP) has a classical solution in [0, $]. 
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Proof. We will give the joint proof for all parts of the theorem, underlying the differences if we 
assume (Ai) or (Ai) 1 . 

Let I be defined by (fT2")) as in Lemma EOl and define operator T by the right-hand side of (JTSJ) , i.e. 

Ty(t):=-l* f Ml)D"<- 2 y(t)dl + D- 2 f(t,y(t))+vot + yo. 



We will consider T acting on Banach spaces C J ([0, a]) endowed with the norms \\y\\i := sup tg [ a ^{\y^ (t)\; < 
j < i}i i = 0, f , 2 and, in particular, Banach spaces 

Q([0,a]):={yeC i ([0,a]);y(0) = 0} 

and 

C5o([0,o]) := {y G C 2 ([0,a]);y(0) = 0,y'(0) = 0} 

endowed with the same norm. We will show that T is a compact operator in quoted spaces and 
that T maps a closed ball into itself. Then we will apply the Schauder fixed point theorem and 
obtain a mild solution. 

(i) T maps C l ([0,a\) into C l ([0,a\) fori = 0,1,2. 

Let y £ C([0, a]), where a comes from (Ai) or (Ai)' . Define 

Mt) ■= / M^D^yWj 



d Mi) 



t 

(t-x) 1 -' r y(x)dxd'y, te[0,a\. (16) 



r(2-7) Jo 

Since y G C([0, a]) and supp^i C [c,d] C [0,2), Jy is continuous. Further, denote 
Ry{t) := (/ * Jy)(i) = / l(x)Jy(t - x)dx, t e [0, a], 



where I is defined as in Lemma 13.31 and therefore is locally integrable. ft follows that Ry is 
continuous. Function Gy defined by 



1-2 



Gy{t) := D- 2 f(t, y{t)) = (t- x)f(x, y(x))dx, t € [0, a] 

Jo 

is continuous by assumption (Ai). Since 

Ty(t):=Ry(t) + Gy(t) + v t + y , te[0,a], (17) 
it follows T : C([0, a]) -» C*([0, a]). 

In what follows we shell need well known fact stated in next lemma. 

Lemma 4.2. (a) Let I € C(K) and suppZ C [0, oo). Then the convolution operator I*, defined 
by y l— > I * y, for y 6 L] oc+ (M), maps continuously C 1 ([0,a]) into Cg([0, a]) and C 2 ([0,a]) into 
C 2 ([0,a]). 

^ Let I G L-J oc+ (R). TTien y I * y, y E Lj oc+ (R) is a continuous mapping from C([0, a]) into 
C ([0, a]), /rom Cq([0, a]) into Cq O ([0, a]) and from C 2 ([0, a]) into C 2 ([0,a]). 



System with distributed fractional derivatives 



9 



Let y G C£([0,a]). We will show that Ty G C£([0,a]). Partial integration in |H} and t/(0) = 
follows 

| ^f^jy x)^C)**r + ,(0) jf 

= JyM(t) te[0,a]. (18) 

Therefore ^</y is continuous and J : Gq([0, a]) — » Gq([0, a]). This, and the fact that Z is locally 
integrable, by Lemma |4T2"1 implies that R maps Gq([0, a]) into Gq([0, a]) and 

^-Ry(t)=l*^Jy(t)=l*JyW(t), te[0,o]. 
at dt 

By (A;) it follows that 

j t Gy(t) = f(x,y(x))dx, te[0,a], 

is continuous and G : Cq ([0, a]) — ► Gq([0, a]). 
In total we have that 

d /"* 

-Ty(t) = I * Jy'(t) + J f(x, y(x))dx + v , t G [0, a] (19) 
is continuous, Ty(0) = 0, and therefore T : Cft([0,a]) -> Cft([0,a]). 

Finally, let y G C$ ([0, a]). We will show that Ty G C? ([0, a}). Since y(0) = and y'(0) = 
d2 r Jy (t) = J„W(t) + y'(0) T -^M^r*, 



rft 2 " v ; " w ' " w y c r(2- 7 ) 



is continuous. Thus, 



+.(o)/ d ^;;W , te[ o, ] 

^Jy(t) = JyM(t), iG[Q,o],*G {0,1,2}, (20) 

Jy(0) = Jy'(0) = 0, so J : Cg ([0,a]) -> C^ ([0,a]). This, with LemmaEjH and the fact that Z is 
locally integrable, implies that R maps Gq O ([0, a]) into Cg O ([0, a]) and 

^Ry(t) = l*-^rJy(t) = l*Jy^(t), te[0,a], ie{0,l,2}. 
Assume now that (A;)' holds. Then 

^Gy(t) = f(t,y(t)), te[0,a] 

is continuous, Gy(0) = ^Gy(0) = and therefore G : C$ o ([0,a]) -> Cg ([0,a]). Thus, T : 
G 2 ([0,a])^G 2 „([0,a]) and 

^Ty(t) = l*Jy"(t) + f(t,y(t)), te[0,a}. (21) 



(ii) Compactness of T. 
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Let us show that J is a compact operator in C([0,a]). Recall, 

N 

\ — 30 ' 



Jy(t) = Jim ^0 1 ( 7 „) J D^- 2 2/ (t)A 7n 



n=l 

,. \- 0i (in) [* y(x) a . rn , 

where 7„ are points of the interval [c, d] in the definition of the Riemann sum defined for the 
integral. Since j n - 2 < d - 2 < 0, for all n E N, operators £H™- 2 : C l ([0,a]) -> ^([O.a]), 
i = 0,1,2, to 6 {1, ...,iV} are operators of fractional integration and therefore compact. J is a 
compact operator as a limit of compact operators in the operator norm topology. By the same 
argument, by ([2"U|) . J : Cq — > Cq and J : Gq — > Gq are compact. 

By Lemma HH1 Z* is a continuous mapping: C([0, a]) — > C([0, a]), Gq([0, a]) — > Cq([0, a]) and 
Cqo([0, a]) — + Cqo([0, a]), respectively, and therefore R is a compact operator as a composition of 
a compact and a continuous operator on C([0, a]), Cg([0,a]) and Gq O ([0, a]), respectively. 
Under the assumption (Ai), G : Cq([0, a]) — ► Gq ([0, a]) is compact. Let M > and Bm be the ball 
in Cft([0,a]) (S A/ := {j/ G Cg([0,o]); ||y||i < M}). The set G[B M ] := {G0;0 G B M } is uniformly 
bounded. Let us prove this. Condition (Ai) implies 



\Gy(t)\ < 

and 



f \t-x\\f(x,y(x))\dx<aM a C, te[0,o], 
Jo 



\j f Gy(t)\< J^\f(x,y(x))\dx<M a C, ie[0,a], 1/6%, 

where C = max te [ a j L* h(x)dx. It is also an equicontinuous family since 
|G7/(t 1 )-G 2 /(t 2 )| < 

/■tl />i2 

</ |(t x - x) - (t 2 - a;)||/(i,j/(i))|da; + / \t 2 - x\\f(x,y(x))\dx 



and 



<M Q G(ii-i 2 )+ sup |t 2 -a:|M a G 

< (h - t 2 )M a C(\ + a) < G'(ti-f 2 ), i!,i a G [0,o] 



I^G&fr) - ^Gfy(*2)| < t \f(x,y(x))\dx 



< M a / \h(x)\dx, h,h G [0,o]. 
Jti 



Since the function i i— > J* |/i(a;)|dic, < G [0, a] is absolutely continuous in [a, 6], it follows that the 
set G[Bm] is equicontinuous. Therefore, Arzela-Ascoli theorem implies that G[Bm] is relatively 
compact and G is, therefore, a compact operator in Gp([0, a]). 

If we use the stronger condition (Ai)' , then G is compact in Gq ([0, a]). We will show that G maps 
a bounded sequence {y n }neN in Gq ([0, a]) to a relatively compact set {Gy n } n ^ in Gq ([0, a]). If 
{?/n}neN is bounded in Gq ([0, a]), then {j/„} n6 n and {i4}neN are equicontinuous and Arzela-Ascoli 
theorem assures that we have convergent subsequences, again denoted by {y n } n eN and {y' n } n eN- 
Denote c„ = Gy n , n G N. Then, for i G {0, 1}, 

d \c n (t) - ^c m (t)\ < f \t - x\^\f(t, y n (t)) ~ f(t,y m (t))\ 



<dt l ny ' dP 

<Di\y n (t)-y m (t)\ a , te[0.a], (22) 
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where Di — max tg [ . a ] J Q \t — x] 1 l \h{x)\dx, i € {0, 1}, and 

l^c„W - ^c m (t)\ < \f(t,y n (t)) - f(t,y m (t))\ 

<D 2 \y n (t)-y m (t)\ a , te[0,a], (23) 

where Di = max (e [ . a ] \h{x)\dx. Since {?/n}neN converges, (|22l) and (|23[) imply that the same is 
true for {c„}„ 6N in C 2 ([0, a]). 

As a sum of compact operators R, G and the mapping: y i— » uq£ + yo 5 is compact and 
T : C[0,a] -> C[0,o]. Moreover, T : Cft([0,a]) -> Cft([0,a]) is compact if (Aj) holds, while 
T : Cg ([0, a]) -> Cg ([0, a]) is compact if (A;)' holds. 

(iii) Determination of 5. 

We will show that there exists (5 > depending on r such that T maps £? r C C([0, 6]) , S r C 
C£ ( [0, 5] ) and B r C Cg ( [0, 8] ) into itself. 
From (fT6|) and ([20]) we derive 

|^Jj/(t)|<M t sup |»W(aO|, t€ {0,1,2}, (24) 
flr xe[o,t] 



where 

^ \Ml)\ 

r(3- 7 )' 

Now consider B r C C{[0,8\) and B r C Cft([0,<5]), respectively. Then JT7J), (HU), and (A/) 
imply that for i = 0, 1, 

||-Tj/W| < |i * ^-Jy(t)\ + \§-Gy(t)\ + + |w,|(l - i) 

t /.t 



< / |l(t-a;)||Jy«(x)|di+ / \(t - x^Wf^y^dx + \vo\t 1 -* + \y \(l - i) 
Jo Jo 

< D t M t sup \y^(t)\ + C t sup \y(t)\ a + {volt 1 " 1 + \y \(l - i), te [0,5], (25) 

where 

D t = I \l{x)\dx, C t = I h{x)(t - xf-t-dx, te[0,5]. 



Since |y^(£)| < r, i = 0, 1 and r > max{|^o|, \yo\}, it is possible to shrink 5 as much as we need 
to reach 

\^Ty(t)\<r, for i — 0,1. 
With this, we determine 5 in cases a) and b). 

For the case c), consider the ball B r C Cq O ([0,S]). Then IpTj). ([24)) and (A/)' imply that the 
estimates for 4riTy(t), i = 0, 1 are valid and 

I^T tf (i)|<|l*^J»(t)| + |/(t,i/(i))| 

< D t M t sup |y< 2 >(t) | + C t \y(t)[ a , te[0,<5], 

with Ct = max tg [o j( s] and -Di and M t as above. Since < r, i = 0, 1, 2 and r > 0, it is 

possible to shrink <5 such that 

\^-Ty(t)\<r, fori = 0,1, 2, 
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so it determines S. 

(iv) Assertions a), b) and c). 

a) Recall that assuming (Ai) and fixing r > max{|y |, |t>o|}j by (i)-(iv), we have obtained the 
existence of 6 — 6(r) > 0, such that T is a compact mapping C([0, S\) — > C([0, (5]) and that T 
maps B r C C([0, <5]) into itself. Thus, according to the Schauder fixed point theorem, T has a 
fixed point in B r C C([0, 5]), i.e. fH} has a solution y S B r C C([0, 5]). Set 



Then (y, z) is a mild solution to (CP). 

b) Again, assuming (Ai) and fixing r, by (i)-(iv), we find 5 = 5(r) > 0, such that T : Cq([0, 5]) — > 
Cq ([0, (5]) is compact and T maps the ball B r C Cq ([0, <5]) into itself. Again, the use of the Schauder 
fixed point theorem assures that T has a fixed point y in B r C Cq([0, <5]). If z is given by (|26p 
then (y, z) is a mild solution to (CP) in [0, 5]. 

If 70 — 7^ > 1, then Lemma [5751 implies that I is absolutely continuous. Thus, I* Jy € AC 2 ([0, a]). 
Also Gy G AC* 2 ([0, <5]). By dTTJ) and the fact that y = Tywe obtain y e AC 2 ([0, a}). 

c) With (A;)' and r > 0, by (i)-(iv), we find <5 = S(r) > such that T is compact and maps 
B r C Cq O ([0, a]) into itself. According to the Schauder theorem, (fH)|) has a solution y in Cq O ([0, 5]). 
If z is given by (f26|) then (y, z) is a mild solution to (CP) in [0, 5]. D 1 y, 7 < 2 is continuous, since 
y G Cq , and z G C([0, <5]) by J2H|). Applying _D 2 to y given by (fTBj) leads to conclusion that (y, z) 
satisfies (CP)\. 



Different conditions on 4>i and / 

Further, we assume restriction on supp^i, which will imply better regularity properties of solu- 
tions. Let 



Then, we have the following result: 

Theorem 4.3. Let yo,vo G K and r > max{|yo|, \vq\}. Let ($2), (Aq) and (Ai) hold. Then there 
exists 5 = S(r) > 0, such that: 

a) (CP) has a non-impact solution in [0,5]. 

b) If (Ai)' holds and yo = 0, then (CP) has a classical solution in [0,(5]. 

Proof. We consider operator T defined by (|17|) and proceed in the same way as in Theorem 14. II 
The proof of this theorem is very similar to the proof of Theorem 14. 1[ therefore we will just give 
the parts of the proof which are different. 

(i) T maps C l ([0,a\) into C J ([0,a]) for i = 1,2. 

Let y G C 1 ([0,a]). The continuity of Ty follows as in Theorem 14. II Let us show the continuity of 
its first derivative -^.Ty. 

An expression for operator Jr J (see ([18])) has an additional summand, which does not vanish (we 
do not have y(0) = 0), but continuity follows from the fact that (^1(7) = for 7 > 1. Therefore, 




(26) 



□ 



($2) 4>i is continuous in [c, d] C [0, 1), ^1(7) = 0, 7 ^ [c, d]; 
4>2 is of the form (fl"T3j) . 



±Ty(t) = ^-(l*Jy)(t)+Gy(t)+vo, 



(27) 
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and thus T : CHIM) -> ^([O.a]). 

Let y G Cq([0, a]) and assume (A;)' instead of (A;). Then, Ty(0) — because y(0) = and the 
first derivative of Ty is given by (|19p . The second derivative given by 

^Ty(t) = J * Jy"(t) + I * y'(0) J J^-t^d^ + f(t, y(t)), t £ [0, a] 

is continuous since y" G C([0, a]), J : C([0, a]) — > C([0, a]), </>i = for 7 > 1 and / is continuous. 
Therefore T : Cg([0, a]) -> ([0, a]). 

(m) Additional regularity properties of T . 

Additional properties of T which assure higher regularity of the solutions are to follow. We will 
show that (Ai) implies the mapping properties 

T : C([0, a}) — > C l ([0, a]) and T : C\[0, a}) -> AC 2 ([0, a]). (28) 

We will use the following lemma. 

Lemma 4.4. aj Let y £ C([0, a]) and a > 1. 7%en J a j/ e C 1 ([0, a]). 
&J Let y G Lj, c ([0,a]) and a > 1. Tften I a y £ AC([0,a]). 

Proof, a) If a = 1 then f t I a y{t) = y(t) and therefore if y G C([0,a]) then I a y G ^([(La]), 
and if y G £^([0, a]) then 7 a y G AC([0,a}). For a > 1 using that I a I^ = I a+0 for a,/3 > 
and D 1 ! 1 = I we have that f t I a y{t) = D 1 I 1 I a - 1 y{t) = I a - l y(t). If y G C([0,o]) then 
J 0-1 !/ G C([0,o]), and again 7 Q y G G 1 ([0,o]). If y G Lj, c ([0,a]) then J"" 1 !/ G L£, c ([0,a]) (since 
/Q_1 : ^locd . a D L loc([°> a D) and therefore J Q y G AC([0, a]). □ 

Let y G C([0, a]). Lemma IOI gives I a y G C 1 ([0, a]) if a > 1, hence 

J : C([0, a]) — > C 1 ([0, a]). (29) 

Since Jy(i) = / </>i(7)D 7 ~ 2 y 0*7, t G [0, a] and d < 1, it follows that operator D 7-2 is in fact / 2 ~ 7 
with 2 - 7 > l; C so L> 7 ~ 2 y G C^fltya]). Since Gy G G^^a]) we have that Ty G G^O,^). 
Let y G C l {[Q,a}). First note that if 7 < 1 then (t 1 " 7 )' = (1 - 7 )t~ 7 G L\ oc {R), hence 

m 1 40if~ idi e ^ c([o ' a]) - (3o) 

Assumption y' G C([0, a]), as well as ([2^) and ([3H|) imply that 

I Jy = Jy' + y(0) | JM-t^dj G AC([0, a]). 

Thus, J : C x ([0,a]) -> AG 2 ([0,a]). Further, Gy G AG 2 ([0,a]) because £?Gy(t) = f(t,y(t)), 
t G [0, a] and for y continuous 1 1— > /(t, y(t)) is a locally integrable function. Therefore, by (|2"T|) we 
have that T : C^a]) -» ,4G 2 ([0,a]). 

(wj T is compact in G l ([0, a]), z = 1,2. 

Let i = 1. By Theorem 14.11 T, given by (jTTJ) (with dT^J)) , maps a bounded sequence in C 1 ([0, o]) 
into a sequence with a convergent subsequence in C 1 ([0, a]). 

For the compactness of the operator T, given by (JTTJ) (with (|27p ). we need additionally to show 
that if {y ra }neN is bounded in C 1 ([0, a]), then 

{Vn(P)J p^%K*)** 1_7 *y}neN 
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has a convergent subsequence. This is true since any bounded sequence {y ra }neN in C 1 ([0, a]) has 
a subsequence, again denoted by {y n }neN, which is convergent in C([0, a}). Hence, {y n (0)}n<EN 
converges as well. 

Similarly, one can prove the compactness of T in Cq ([0, a]). 
(iv) Determination of 8. 

Take P r C ^([0,(5]) and P r C Cq([0,S]), respectively, and 6 will be defined later. Then, for 
^G{0,1}, 

\l * 2/ W (0) | f^Z^ 1 " 7 ^! < |j/ (i) (0)|M t ' < M t '||y|| m , t G [0, <5], 



where 



I// : / /lr),/r / J|M_ t i-7 d7i teM ]. 



Using this, (|2"T|) and estimate (JUJ), similarly as in the proof of Theorem 14. 1\ one can estimate 
4pTy(t), for i = 0, 1, if (A;) is assumed, and for i = 0,1,2, if (Ai)' is assumed. Therefore, it is 
possible to find <5 > such that y G B r C C 1 ([0, <5]) and y £ B r C Cq ([0, <$]), respectively, implies 
that 

h^-ry(0|<r, i = 0,l, respectively i = 0,1,2. 
(v,) Assertions a) and b). 

a) The Schauder fixed point theorem implies the existence of a fixed point y for T in C 1 ([0, <5]), if 
(A;) is assumed. Let z be given by (|2U1) . Then (y, z) is mild solution to (CP). By (|28p we obtain 
that y G AC 2 ([0,<5]). This, with ef < 1, implies 

d 

M7)D y y(t)d 7 G AC([0,<5]). 

Since £ is locally integrable, we have that z G AC([0, <5]). Therefore, the mild solution (y, z) satisfies 
the first equation in (CP) in Ll oc ([0,5]), so this is also a non-impact solution 

b) If we assume (Ai)' , then T has a fixed point in Cq([0,$\), (CP) has a mild solution (y, z) 
in [0,5], z is given by and y G Cq([0, <5]). So (CP) 2 holds. Further, rf < 1 implies that 
f d Ml) D1 y(t)dl e ^([0,(5]) and therefore, 2 G C x ([0,<$]) and (CP)i is satisfied for all t G 
{0,6}. " □ 

In this particular case, one can also impose a condition on /, stronger then (Ai) and (Ai)' , and 
obtain classical solution to (CP) in [0,(5], which satisfies (CP) 2 for all i G [0,(5]. The condition 
reads: 

(Ai)" f(t,u), (t,u) G [0, oo ) x M, is continuous f(t,0) = 0,t G [0,oo), and 
there exist a > 0, ft € C([0, a]) with ft, > and ft(0, 0) = 0, such that 
\f(t,u) - f(s,v)\ < h(t,s)(\t - s\ + \u- v\), u,v G R,t G [0,a]. 

The corresponding theorem is the following. 

Theorem 4.5. Let vq Gl, r > |uo| and yo — 0. Let ($2), (A)) a ^ (A)" hold. Then there exists 
(5 = S(r) > 0, sitcft that there exists the classical solution to (CP) in [0,(5], which satisfies equation 
(CP) 2 point-wisely for t G [0, <5]. 



Proof. First note that all assumptions of Theorem 14.31 arc satisfied. Hence there is a classical 
solution to (CP) in [0,(5] obtained in five steps in previous proof. The assumption (A{)" gives 



System with distributed fractional derivatives 



15 



stronger results. Note that if y G AC([0, a]) and (Ai)" holds, then t t— > f(t,y(t)), t G [0, a] is 
absolutely continuous. Indeed, 

\f(t,y(t)) - f(s,y(s))\ < h(t,s)(\t -s\ + \y(t) - y(s)\) 

implies absolutely continuity since y belongs to AC([0, a}). 

Further note that in step (ii) of the proof of the Theorem 14.31 in addition to l|28p we have 

T : C 2 ([0, a]) — > AC 3 ([0, a}). (31) 
To see this, let y G C$([0, a]). Then dHJ) and imply that 

Since £ i— » /(i, ?/(£)), i G [0, a] is absolutely continuous, we have that Gy G j4C 3 ([0, a}) and therefore 
(l3"Tj) . It follows that the fixed point y G Cq ([0, a]) obtained in part b) of the fifth step of the proof 
of the Theorem 14.31 is then an element of AC 3 ([0, 5}). Then z given by is an element in 
AC 2 ([0, S]), which yields that both equations in (CP) are satisfied for all t G [0, 6}. □ 

Remark 4.6. a) Note that in Theorem 14.11 and Theorem 14.31 we could have assumed that [c, d] 
was a subset of the interval (— oo, 2) and (— oo, 1), respectively, in which case we would have the 
same results. Indeed, with respective assumptions, we change the bounds of the integral in (fTS)) . 
the definition of operator J, but not its properties: jhjJy, % = 0, 1, 2 remain continuous, J remains 
compact in respective spaces, and estimates in (|23|) hold. 

b) With assumptions d < 2 and d < 1, respectively we could also have considered a constitutive 
equation of the form 

and we would obtained the same results as in Theorem 14.11 and in Theorem 14.31 respectively. 

c) It also make a sense to consider constitutive equation for d < 0. Then on its right hand side 
only fractional integrals of y appears. Such case is covered with condition that follows. 



($3) <pi is continuous function in [c, d] C (— 00, 0), ^1(7) = 0, 7 ^ [c, d]; 
<f>2 is of the form (|13[). 

Also, the theorem similar to Theorems 14.11 and 14.31 holds. 

Theorem 4.7. Let yo,VQ G M and r > max{|yo|, \ v o\}- Let ($3), (Aq) and (Ai)' hold. Then 
there exists S = 5(r) > such that (CP) has a classical solution in [0,8]. Moreover, the classical 
solution satisfies (CP)2 for all t G [0,5]. 

Finally, we are interested in cases when both (j>\ and (f>2 are linear combinations of translations of 
delta distributions: 

m 

($4) M-) =£&i*(- e [0,2), 0o > 0i > m ,j G {0,...,m} 

4>2 is of the form (fl~3l) . 
The constitutive equation (CP)2 becomes 

k rn 

Y,a i D^z(t)=J2bjD l3j y(t), <>0, 

i=0 j=0 

and the theorem similar to previous holds. 
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Theorem 4.8. Let yo,vo £ R and r > max{|y |, |^o|}- Let ($4), (A)) and (A;) Wd. Tften f/iere 
exists 5 — S(r) > suc/i i/ia£; 

a) (CP) has a mild solution in [0,(5]. 

b) (CP) has a non-impact solution in [0,(5] if one of the following conditions hold 

(i) yo = and 70 - 7& > 1 or (m) /3 < 1. 

If (Ai)' holds, then (CP) has a classical solution in [0,(5] if one of the following conditions 
hold 

(i) yo = vq = or (ii) fio < 1 and yo = or (Hi) /3q < 0. 
The proof similar to previous ones is skipped. 

5 Continuous 02 

In this section we impose different condition to <p2- 

($5) : • (/>i is continuous function in [c, d] C [0, 2), (^1(7) = 0, 7 ^ [c, d] or 

m 

01 (0 = - A')' A e [°' 2 )>Ai > ft > Pm,3 G {0,...,m} 

3=0 

• 02 G C 3 ([0, 1]), (f> 2 = out of [0, 1], </» 2 (1) / and 
either 0(0) ^ or (f> 2 (7) — P7 9 , p > 0, g > 0. 

Assumptions on </>2 given in ($5) will imply (as it is shown in |18j ) the existence of a solution to 
distributed order differential equation J Q 4'2(l)D 1 y — g, i.e. the existence of a locally integrable 
function C^ 1 ( ri , } . ^ , ). 

Theorem 5.1. Let yo,Vo el, r> max{|?/o|, |uo ] } an d f satisfy (Ai). Let <f>\ and (f>2 satisfy ($5). 
Then there exists 5 — S(r) > such that: 

a) There exists a mild solution of (CP) in [0,(5]. 

b) There exists a non-impact solution of (CP) in [0,(5] i/supp0i C [c, 1). 

c) There exists a classical solution of (CP) in [0,5] if (A{)' and one of the following conditions 
hold: 

(i) yo = vq = or (ii) supp</>i C [c, 1) and yo — or (Hi) supp</>i C [c, 0). 



Proof. The solution to equation 



1 

02(7)-D 7 zdj = g, in S?\. 



is given by 

z = X * 9, where x = C~ 1 (— r —— -), 

Jo <t>2 (7) « 7 07 

provided that the inversion exist. It is proved in O Proposition 3.1], that conditions on <f>2 given 
in ($5), imply that x exists and, moreover, that x G C°°((0, 00)) nL 1 1 oc ([0, 00]). Therefore, (CP) 2 

1 ,2 

fairfD^zdj = / (j) 1 (j)D 1 yd'f 
Jo Jo 

has a solution in 

z = x* (j)i(i)D 1 yd-y. 
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The substitution of z in (CP)\ and the integrations give 

V = -X* [ (l>{lW- 2 yd 1 + D- 2 f{- 1 y{-)) + v Q t + y , in 



Again, we consider operator T which is given by (|17p . where R is now given by Ry := x* Jy and 
J by (|16p . Since x 1S locally integrable, the same properties hold for J, R and T, as in Theorems 
14.11 and others. Thus, it follows that there exist 5 = S(r) and a mild solution (y, z) to (CP) in 
[0, S}. Also, one can prove the properties of y and z quoted in a), b), and c) in the same way as in 
Theorems 14.11 and others. 

□ 

To close the section we give the example that follows arise from application. Consider the system 
(CP) with </>i(7) = 6 7 and ^2(7) = a 7 7 7 £ [0, 2), where a and b are positive constants with b > a. 
The latter condition is consequence of the Second Law of thermodynamics. 



(CP)e 



D 2 y(t) + z(t) = f(t,y(t)), t>0, 
WD"<y(t)d'y = f* a^D^z(t)d-y, t > 

y(0) = y o ; y( 1 )(0) = v , 

One can apply Theorem 15.11 to prove the existence of a mild and a classical solution for above 
system. 

Theorem 5.2. Let yo,vo <E R and r > max{|yo|, l^ol}- Let (Ai) hold. Then problem (CP) ex 

(i) has a non-impact solution. 

(ii) has a classical solution (y,z) if (A{)' holds and if yo = 0. 

Proof. Since X ,0 2 £ C 3 ([0, 1]), (f> 2 (l) = b ^ and <f> 2 (0) = 1^0 conditions of TheoremO 
(part b) case (ii) and part c) case (ii)) are satisfied. Theorem 15. II implies the assertions. □ 
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